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ABSTRACT 



Aims. We prepare real-life Cosmic Microwave Background (CMB) lensing extraction with the forthcoming Planck satellite data, 
by studying two systematic effects related to the foregrounds contamination: the impact of foreground residuals after a component 
separation on the lensed CMB map, and of removing a large contaminated region of the sky. 

Methods. We first use the Generalized Morphological Component Analysis (GMCA) method to perform a component separation 
within a simplified framework which allows a high statistics Monte-Carlo study. For the second systematic, we apply a realistic mask 
on the temperature maps and then, restore them using a recent inpainting technique on the sphere. We investigate the reconstruction 
of the CMB lensing from the resultant maps using a quadratic estimator in the flat sky limit and on the full sphere. 
Results. We find that the foreground residuals from the GMCA method does not alter significantly the lensed signal, nor does the 
mask corrected with the inpainting method, even in the presence of point sources residuals. 

Key words. Cosmic microwave background - Gravitational lensing - Large-scale structure of Universe - Methods: statistical 



1. Introduction 

Cosmic Microwave Background (CMB) temperature 
anisotropies and polarisation measurements have been one 
of the key cosmological probes to establish the current cosmo- 
logical constant A and Cold Dark Matter (ACDM) paradigm. 
Reaching the most precise measurement of these observables is 
the main scientific goal of the forthcoming or on-going CMB 
experiments - such as the European Spacial Agency satellite 
PlanckQ which has been successfully launched on the 14th of 
May 2009 and has currently begun collecting data. 

Planck is designed to deliver full-sky coverage, low-level 
noise, high resolution temperature and p olarisations maps (see 
iTaubetfeOOfjtlThe Planck Consortial2005T) . With such high qual- 
ity observations, it becomes doable to extract from the CMB 
maps cosmological informations beyond the angular power 
spectra (two-points correlations, hereafter A PS), by exploitin g 
the measurable non-Gaussia nities (see e. g. iKomatsul d2002l) ; 
iThe Planck Consortia! d2005l) ). 

The weak gravitational lensing is one of the sources of 
non- Gaussianity affecti n g the CMB after the recombination 
(see iLewis & ChallinoJ I2OO6I) for a review). The CMB pho- 
tons are weakly deflected by the gravitational potential of 
the intervening large-scale structures (LSS), whi ch perturbs 
the Gaussian statistic of the CMB anisotropies dBernardeaul 



1 19971: IZaldarria"gal 120001) . Conversely, it becomes possible to 
reconstruct the underlying gravitational potential by exploit- 
ing the higher-order co rrelations induced by the weak lens- 
ing in the CMB maps dBernardeaul 1 1997t iGuzik et al J 120001: 
iTakada & Futamasdl200UlHdl2001btlHirata & Seliakll2003al) . 



The relevance of the CMB lensing reconstruction for the 
Cosmology is twofold. First, for the sake of measuring the pri- 
mordial B-mode of polarisation predicted by the inflationar y 
models dKamionkowski et al.ll997tlSeliak & Zaldarriagal 19971) . 
the CMB lensing is a major contaminant. It induces a sec- 
ondary B-mode p olarisation signal in perturbi ng the E-mode po- 
larisation pattern dZaldarriaga & Seliakll 1998b . A lensing recon- 
struction allowing the delensing of the C MB maps is required 
to recover the primor dial B mode signal dKnox & Sond 120021: 
ISeliak & Hiratall2004l) . However, CMB lensing is also a power- 
ful cosmological probe of the matter distribution integrated from 
the last scattering surface to us. In a near future, it is an unique 
opportunity to probe the full-sky LSS distribution, with a maxi- 
mum efficiency at redshift around 3, where structures still expe - 
rience a well described linear growth dLewis & Challinoj2006l) . 
A lensing reconstruction would largely improve the sensitivity 
of the CMB experiments to the cosmological parameter affect- 
ing t he growth of the LSS, such as neutrino mass or dark en- 



- g ri 

ergy dHull2002tlKaphnghat et alJl2003l:lLesgourgues et alj|2006t 
iPerotto et al.ll2006l) . 
Although well-known theoretically (Bl anchard & Schneidej 



1 http : //www . r ssd . esa . int/index . php?pro j ect=PLANCK 
&page= index 



1987 1), the CMB lens ing has never bee n directly measured. 
ISmith et all (120071) and lHirata et al.l d2008l) have found evidence 
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for a detection of the CMB lensing in the WMAP data by corre- 
lating with several other LSS probes (Luminous Red Galaxies, 
Quasars and radio sources) at 3.4<x and 2.5<x level respec- 
tively. This situation is expected to change with the forthcom- 
ing Planck data. Planck will be the first CMB experiment al- 
lowing to measure the underlying gravitational potential without 
requiring any external data. However, even with the never be- 
fore met quality of the Planck data, CMB lensing reconstruc- 
tion will be challenging. CMB lensing is a very subtle sec- 
ondary effect, affecting the smaller angular scale at the limit 
of the Planck resolution in a correlated way over several de- 
grees on the sky. As already quoted, CMB lensing reconstruc- 
tion is based on the induced non-Gaussianities in the CMB 
maps, in the form of mode coupling. Consequently, any pro- 
cess resulting in coupling different Fourier moments is a chal- 
lenging sy stematic to d e al wit h in order to retrieve the lensing 
signal (see ISu & Yad av (2009) for a recent study of the impact 
of instrumental systematics on the CMB lensing reconstruction 
bias). Astrophysical components and other secondary effects 
might also be a source of non-Gaussianity. These components 
include : Thermal and Kinetic Sunyaev-Zel'dovich effects (thSZ 
and kSZ), due to the scatte ring of CMB radiation by el ectrons 
within the galaxy clusters (ISunvaev & Zeldovichl 119701) : fore- 
ground emissions, such as synchrotron, Bremsstrahlung and dust 
diffuse galactic emission as well as extragalactic point sources. 
All these components may give a sizabl e contribution to the level 
of non-Gaussianities in the CMB maps (lAghanim & Fornil 19991: 
Argiieso et al l 120031: lAmblard et al.N20o4 ; iRiquelme & Spergell 



20071: iBabich & Pierpaolill2008l) . 



The impact of most of the aforementioned effects on the 
C MB lensing analysi s with WMAP data has been investigated 
bv lHirata et al.l (12.0081) . The fact that they found a negligible con- 
tamination level is encouraging. However, such a result could 
change when one considers the high er resolution, better s ensi- 
tivity maps provided by Planck. In iBarreiro et al.l (120061) . the 
component separation impact on non-Gaussianity was studied 
in the framework of the PLANCK project, but no lensing recon- 
struction was performed. Hence, the impact of these foreground 
residuals on the CMB lensing reconstruction is still to be stud- 
ied. 

The overall purpose of the present study is to give an in- 
sight to the issues we should deal with before undertaking any 
complete study of the CMB lensing retrieval with Planck: what 
is the impact of the foreground residuals on the CMB lensing 
reconstruction? Will it still be possible to reconstruct the CMB 
lensing after a component separation process or will such a pro- 
cess alter the temperature map statistics? How should we deal 
with the masking issue? Beyond the detection of the CMB lens- 
ing signal, we tackled the reconstruction of the underlying pro- 
jected potential APS. We investigated two issues, the impact of 
a component separation algorithm on the lensing reconstruction 
and the impact of a masked temperature map restoration before 
applying a deflection estimator. 

Sect. |2]briefly review the CMB lensing effect and the recon- 
struction method. We present in Sect. [3] an analysis of the im- 
pact of one component separation technique, named Generalized 
Morp hological Component Analysis (GMCA) dBobin et al.l 
2008), which is one of th e different method s investigated by 
the PLANCK consortium dLeach et al.l 120081) . In Sect. |4] we 
show how a recen t gap filling method (i.e. inpainting process) 
dAbrial et al.ll2 008) may solve the masking problem, which may 
be one of a the major issue for the CMB lensing retrieval since 
it introduces some misleading correlations between different an- 
gular scales in the maps. 



2. CMB lensing 

In this section, we briefly review the CMB lensing effect and 
the reconstruction method. We introduce the notations used 
throughout this paper. 

CMB photons geodesic is weakly deflected by the grav- 
itational potential from the last scattering surface to us. 
Observationally, this effect results in a remapping of the 
C MB temperature anisotropics T = A0/0cmb, according 
to lBlanchard & Schneider] d 1987b : 



T(h) = T(h + d(h)). 



(1) 



In words, the lensed temperature T in a given direction of the 
sky h is the temperature T one would have seen in the neighbor- 
ing direction h + d(h) in the absence of any intervening mass. 
The deflection angle, d(h), is the gradient of the line-of-sight 
projection of the gravitational potential, d{h) = V0(«), where 
<p can be calculated within the Born approximation, as the in- 
tegral alo ng the line-of-sight of the tridimensional gravitational 
potential dChallinor & Lewisll2005l) . 

The CMB lensing probes the intervening mass in a broad 
range of redshifts, from z* = 1090, at the last scattering sur- 
face, to z = 0, with a maximum efficiency at z ~ 3. At such 
a high redshift, the LSS responsible for the CMB lensing (with 
typical scale of 300/j -1 M pc) still experience a linear regime of 
growth. As a result, the projected potential <p can be assumed to 
be a Gaussian random field; The consequences of the non-linear 
correction s to c/> are shown to be w eak on the CMB lensed ob- 
servables (Ch allinor & Lewisll2005l) . Thus, this hypothesis holds 
very well as long as the CMB lensing study doesn't aim at mea- 
suring a correlation with other LSS probes at lower redshifts. 

Besides, in the standard ACDM model, the deflection angles 
have a rms of - 2.7 arcmin and can be correlated over several 
degrees on the sky. The typical scales of the lensing effects are 
small enough for a convenient analysis within the fiat sky ap- 
proximation. The projected potential may be decomposed on a 
Fourier basis <p{k). and its statistics is completely defined by: 



< 0(Jfci)0(*2) >= (2;r) 2 5(A: 1 + Jfca) C* 



(2) 



where C ki is the full-sky projected potential APS taken at a mul- 
tipole / = \k\ | and it is related to the deflection APS by : 



(3) 



The lensed CMB temperature APS can be derived from th e 
Fourier transform of Eq. Q (e. g. as in Ok amoto & Hul (120031) ). 
Lensing effect slightly modifies the APS of the CMB tempera- 
ture, weakly smoothing the power at all angular scale at the ben- 
efit of the smaller angular scales. Deeply in the damping tail, at 
multipole / > 3000, lensing contribution even dominates over the 
pure CMB one. However, the main observational consequences 
of the CMB lensing effect lies beyond the APS. The remapping 
induces non-Gaussianities in the CMB temperature field, in the 
form of some correlations between different angular scales. 

Consequently, the two-point correlation function of the 
lensed temperature modes, calculated at the first order in 0, 



2 A priori, the remapping function should depend not only on a con- 
vergence field but also on a rotation field, so that the deflection an- 
gle is not purely gradie nt but have a rotational contribution. However, 
iHirata & Seliakl j2003bl) have shown that rotation field effect would be 
negligible for the next generation of CMB experiments. 
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writes dOkamoto & Hull2003l) : 

<f(*: 1 )f(A: 2 )>cMB = (2nfS(k x + k 2 )C^ 

+ frrikukoJ^m + Oiip 2 ) 



(4) 



where L = k\ + k 2 , and the CMB subscript denotes an ensemble 
average over different realisations of the CMB but a fixed inte- 
grated potential field. The weighting function fjj depends on the 
primordial temperature APS, such as: 



fn(k u k 2 ) - L k\ C™ + Lk 2 C" 



(5) 



Similarly, one can calculate the four-poin t correlation functio n 
of the CMB temperature field - such as in iKesden et al.l d2003l) . 
One finds that the trispectrum of the lensed temperature field - 
or equivalently, the connected part of its four-point correlation 
function - is non-null even if the underlying (unlensed) temper- 
ature field is purely Gaussian. 

For the sake of reconstructing the integrated gravitational po- 
tential field from a lensed CMB map, two methods exist in the 
literature. Firs t , the quadratic estimator appro a ch has been de- 
veloped in lHul d2001bl) : liiu & Okamotol d2002l) ; lOkamoto & Hul 
d2003l) . Then a maximum-likelihood estimator method has been 
derived in lHirata & Seliakl (|2003a b). The latter approach may 
increase the capabilities of the highest sensitivity highest reso- 
lution CMB projects in reconstructing the integrated potential 
while the former method is still very close to the optimality 
for current built experiment such as Planck. Thus we adopt the 
quadratic estimator throughout this work. 

In the flat sky appro ximation, the estimate d potential map 
takes the following form dOkamoto & Hull2003l) : 



TT (L) = 



Att(L) 
L? 



J (2tt) : 



T(k l )T(k 2 )F TT (k l ,k 2 ), 



(6) 



where the Fourier modes T(k) refer to the observed temperature 
modes, affected by both the CMB lensing and the instrumen- 
tal noise of the CMB experiment concerned. More precisely, the 
temperature map is assumed to be contaminated by an additional 
white Gaussian noise and deconvolved from a beam function as- 
sumed to be Gaussian, so that its APS reads: 



(7) 



where NT t is the instrumental noise APS, modeled in this anal- 



ysis as: 



< 2 ^fwhm 



81n2 



(8) 



where #f w h m is the full-width at half maximum (FWHM) of the 
beam function and <x T , the root mean square of the noise per 
resolution elements. 

Besides, the normalisation function is calculated so that 
(pjj(L) is an unbiased estimator of the integrated potential field: 



Att(L) = L 



d 2 ki 



frjikuk^Frrikuko) 



(9) 



Then, the weighting function Fjj is adjusted to minimize the 
dominant contribution to the estimator variance, i. e. the uncon- 
nected par t of the quantity TT (L )4> TT (L')) - {2nf 5{L + L')Cf. 
Derived in lOkamoto & Hul d2003l) . the calculation leads to: 



F TT (ki,k 2 ) = 



fn(k u k 2 ) 



(10) 



where Cj T = Cj 1 +N^ T is the observed temperature power APS 
as defined in Eq. Q. 

Finally, the covariance of the integrated potential field esti- 
mator provides us with a four-point estimator of the integrated 
potential APS. When expanding the lensed CMB temperature 
modes at second order in 0, the <pjj estimator covariance reads: 



*tt(L)0tt(L')> = (2k) 2 6(L + L')C^ t JL), 



(11) 



where the estimated potential APS, C^ TT (L), taking into ac- 
count all sources of variances, both projected potential and CMB 
cosmic variance, instrumental noise and confusion noise from 
other potential Fourier modes, writes: 



c£ iTT (L) = Cf + N^(L) + N»(L) + N»(L). 



(12) 



Here, we have distinguished three different noise contributions 
to the integrated potential estimator variance. The dominant 
noise contribution, N^(L) = Ajj(L), depends only on the 
unlensed and the observed temperature APS. It represents the 
Gaussian contribution to the potential APS estimator, in a sense 
that it is the variance one would obtain by replacing the lensed 
temperature map in Eq. d5]l, by a map with the same APS but 
Gaussian statistics. In addition, the potential APS estimator suf- 
fers from sub-dominant non-Gaussian noi se contributions. Th e 
first, quoted N**(L) has been calculated bv lKesden etaT] d2003l) . 

The second non Gaussian noise te rm Nt ^jL) is quadratic in C L . 
First calculated in lHanson et al.l d2009l) . this term is shown to 
contribute even more than the first-order one at low multipoles 
(L < 200). 

These two non Gaussian noise terms arise from the trispec- 
trum part (or so-called connected part) of the four-point lensed 
temperature correlator hidden in the integrated potential field es- 
timator covariance. It can be interpreted as the confusion noise 
coming from other integrated potential modes. Since it depends 
on the integrated potential APS, which has to be estimated, an 
iterative estimation scheme would be required for taking it into 
account. However, our study based on simulated data allows us 
to calculate these terms from the fiducial potential APS and then 
subtract it from the estimator variance. 

3. Effect of foreground removal: a Monte-Carlo 
analysis 

Up to now, no analysis has been performed to assess the effect of 
a component separation process on the CMB lensing extraction. 
The question we propose to address here is whether the lensing 
signal is preserved in the CMB map output by the component 
separation process. In order to get a first insight, we use a Monte- 
Carlo approach within the flat sky approximation. 

3. 1 . Idealized Planck sky model 

We create a simulation pipeline to generate some idealized 
synthetic patches of the sky for the Planck experiment. Our 
sky model is a linear uncorrelated mixture of the lensed CMB 
temperature and astrophysical components, which includes the 
Sunyaev-Zel'dovich effect, the thermal emission of the interstel- 
lar dust and the unresolved infrared point sources emission. In 
modeling these three components, we ensure to catch the dom- 
inant foreground emission features at the Planck-HFI frequen- 
cies. Then we add the nominal effects of the Planck-HFI instru- 
ment, modeled as a purely Gaussian shaped beam and a spa- 
tially uniform white Gaussian noise. Each hypothesis we adopt 
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is a crude model of the astrophysical contaminant and system- 
atic effects that pollute the Planck data, and is intended to be 
a demonstration model for a study devoted to the impact of the 
component separation algorithms on the CMB lensing retrieval. 

We generate four sets of 300 Planck-HFI synthetic patches 
of the sky, with instrumental noise and, when needed, with fore- 
ground emissions: 

- Set I contains lensed CMB temperature maps generated from 
an unique fixed projected potential realisation and with the 
instrumental effects (beam and white noise); 

- Set I-fg is built from set I. In addition, a fixed realisation of 
dust and SZ is added to each set I map; 

- Set II is a set of lensed CMB temperature maps generated 
from 300 random realisations of the lenses distribution plus 
the instrumental effects; 

- Set Il-fg is built from set II. Each map of set II is superim- 
posed with randomly picked dust and SZ maps. 

Note that the point sources emission will be included after- 
ward in our simulation pipeline, through a direct estimation of 
the point sources residuals after component separation as de- 
scribed in Sect. 13.21 Sets I and I-fg will serve at studying the 
projected potential field reconstruction, whereas sets II and Il-fg 
will be used in the projected potential angular power spectrum 
(APS) estimates analysis. Hereafter, our method and its assumed 
hypothesis are detailed. 

3.1 .1 . Lensed CMB temperature map 

Once we have assumed the Gaussianity of the integrated poten- 
tial field, the lensed CMB temperature simulation principle is 
straightforward as a direct application of the remapping Eq. ([TJ. 
We start from the APS of both the temperature and the projected 
potential field as well as the cross-APS reflecting the correla- 
tion between the CMB temperature and the gravitational poten- 
tial fields due to the Integrated Sachs-Wolfe (ISW) effect. Then 
we generate two Gaussian fields directly in the Fourier space, so 
that 



(13) 



p TT 



(O 2 



1*1 



1*1 - 



t (k) and G^\jk) are two independent realisations of 



T(k) = 
<P(k) = 
where G 

a Gaussian field of zero mean and unit variance. Because of the 
typical scales of the deflection field - deflection angles are of 
the order of 2 or 3 arcmin (depending on the fiducial cosmo- 
logical model) but correlated over several degrees on the sky 
- the generated maps should be both high resolution and ex- 
tended over not too small sky area. We choose to produce some 
12.5 x 12.5 square degrees maps of 2.5 arcmin of resolution, as 
a good trade-off between the quality of the simulated maps and 
the time needed to the generation and the analysis of these maps. 

From CMB temperature and projected potential in the 
Fourier space, we calculate both the temperature and the deflec- 
tion angles in the real space. The last step consists in performing 
the remapping of the primordial temperature map according to 
the deflection angles. Here is the technical point. Starting from 
a regular sample of a field (the underlying unlensed map), we 
have to extract an irregular sample of the same field (the lensed 
map) - the new directions where to sample from are given by 




Fig. 1. CMB temperature APS. Red/black (respectively 
green/grey) line is the lensed (respectively unlensed) tem- 
perature APS calc u lated w ith the public Bo lt zmann code 
camb dLewis et~a D l2000t IChallinor & Lewis! 120051) . The 
blue/black data-points are the mean of the binned power 
spectrum reconstructed on 500 simulated lensed temperature 
maps. The error-bars are given by the variance of the 500 APS 
estimates. 



the previous one shifted by the deflection angles. Thus, this is a 
well-documented interpolation issue, the difficulty lying in the 
fact that the scale of the interpolation scheme is the same than 
the typical scale of the physical process under interest. We have 
to take particular care in the interpolation algorithm to avoid 
creating some spurious lensing signal or introducing additional 
non-Ga ussianities. We find that a par ametric cubic interpolation 
scheme dPark & Schowengerdtlll 983*1) apply reasonably well. In 
addition, to avoid any loss of power due to the interpolation, 
we overpixellize twice the underlying unlensed temperature and 
deflection field. The first test we perform to control the qual- 
ity of the simulation is to compare the Monte-Carlo estimate of 
the APS over 500 simulations of the lensed maps with the an- 
alytical calculatio n of the lensed APS provided by the cami 3 
Boltzmann code dLewis et al. 1 120001: IChallinor & Lewis! I2005TT 
As shown in Fig. [1] the APS of our simulated lensed maps is 
consistent with the theoretical one up to multipole 4000 - which 
is largely enough to study the CMB lensing with Planck. 



3.1.2. Astrophysical components 

In any CMB experiment the temperature signal is mixed with 
foreground contributions of astrophysical origin - among them 
we can separate the diffuse galactic emission (thermal and ro- 
tational dust, synchrotron, Bremsstrahlung (free-free) radiation) 
from the extragalactic components (point sources, thermal and 
kinetic Sunyaev-Zel'dovich effects). As discussed in the intro- 
duction, each of these components could potentially, if ineffi- 
ciently removed, degrade our capability to reconstruct the CMB 
lensing. Here, to complete our demonstration sky model, we 
choose to simulate the dominant astrophysical foregrounds at 
the Planck-HFI frequencies, namely the thermal emission of the 
galactic dust, the thermal SZ effect and the unresolved infra- 
red point sources. Thermal dust simulations are obtained from 
an interpolation of the 100;um IRAS data, in the sky region lo- 
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cated around a = 204° a nd 6 - 11°, at the r elevan t CMB fre- 
quencies as described in iDelabrouille et alj d2003l) . Note that 
several treatments have been applied on these maps - point 
sources removal, destriping, inpainting in the Fourier space with 
constraint realisations - which may induce an amount of addi- 
tional non-Gaussianities. The SZ emission on the sky patches 
can be randomly selected in a set of 1500 realisations produced 
with a semi-analytical simu lation tool provided in the litera- 
ture ( IDelabrouille et al.l2002l) . Note that SZ emission is assumed 
here not to correlate with the CMB lensed signal. As for the es- 
timates of the unresolved point sources residuals, we choose to 
take advantage of the refined full-sky simulations of the infra- 
red point sources emission in each Planck-HFI frequencies, 
provided by the Planck Component Separation Working Group 
(WG2). In these simulations, the source counts are drawn from 
the IRAS catalog, and their spectral energy distributions are 
modeled following lSerieant & Har rison (2005). In addition, they 
involve several refinements, such as the filling of the IRAS mask 
by synthetic dat a, the additional sim ulation of fainter sources 
accordin g to the | Granato et alj d2004l) model, and their cluster- 
ing (see iLeach et all d2008l) ). Note that the radio-galaxies, an- 
other population of extragalactic sources, can be safely neglected 
here, as they lead to a sub-dominant emission compared to the 
infrared-galaxies one, at the Planck-HFI frequencies. 

3.1.3. Planck- I ike noise 

Finally, we simulate the effects of the Planck High Frequency 
Instr ument (HFI) according to their nominal characteris- 
tics dThe Planck Consortial [2005), which are summarized in 
Table Q] At each frequency channel, the component mixture is 
convoluted by a Gaussian beam with the corresponding FWHM 
size. Then a spatially uniform white noise following a Gaussian 
statistic is added. Finally, the resulting maps can be deconvolved 
from the beam transfer function, resulting in an exponential in- 
crease of the noise at the scales corresponding to the beam size. 
Because smaller angular scale carry the larger amount of lens- 
ing information, the higher the angular resolution is, the bet- 
ter the lensing reconstruction can be. Our tests show that in 
the ideal case, the lensing reconstruction on Planck-HFI syn- 
thetic maps is insensitive to the addition or the removal of the 
100 GHz frequency channel information, whose beam function 
is roughly twice larger than the beam in the higher frequency 
channels. Worst, when running the full Monte-Carlo chain, turn- 
ing on the foreground emissions and the component separation 
process, adding the lower frequency channel results in increasing 
the confusion noise of the lensing reconstruction. That's why, we 
exclude the 100 GHz frequency channel of our analysis. To sum- 
marize, our Planck sky model reads: 

T obs =y[s + FT 1 ®n, (14) 

where T obs is the set of five individual frequency channel maps, 
y{ is the mixing matrix calculated from the frequency depen- 
dence of the three signal contributions, s is the set of CMB, dust 
and SZ maps at a pivot frequency, B~' is the set of Gaussian 
beam inverse transfer functions and n, the set of white noise 
maps in the five frequency channels we have selected. Note that 
we assume a perfect beam deconvolution process. 

3.2. Component separation using GMCA 

For most of the cosmological analysis of the CMB data - and 
for the CMB lensing extraction in particular - the cosmological 



5 

Table 1. Instrumental characteristics of Planck-HFP 



channel (GHz) 


fwhm (arcmin) 


cr T (/jK-.sr ) 


100 


9.5 


6.8 


143 


7.1 


6.0 


217 


5. 


13.1 


353 


5. 


40.1 


545 


5. 


401 


857 


5. 


18 300. 



a (see lThe Planck Consortia 2005) 



signal has to be carefully disentangled from the other sources of 
emission that contribute to the observed temperature map. The 
component separation is a part of the signal processing dedi- 
cated to discriminate between the different contributions of the 
final maps. Briefly, the gist of any component separation tech- 
nique consists in taking profit of the difference in the frequency 
behavior and the spatial structures (i.e. morphology) that distin- 
guish these different components. From a set of frequency chan- 
nel maps, a typical component separation algorithm provides a 
unique map of the CMB temperature with the instrumental noise 
and a foreground emission residual. In general, the lower the 
foreground residual rms level the better the separation algorithm. 
However, this simple rule is not necessarily true for CMB lens- 
ing reconstruction. In this case, preserving the statistical proper- 
ties of the underlying CMB temperature map is critical. 

In the Planck consortium, the component separation is a crit- 
ical issue, involving a whole Working Group (WG2) devoted 
to provide several algorithms for separati ng CMB from fore- 
grounds and to compare their merits (see ILeach et"aT] d2008[) 
for a recent comparison of the current proposed methods). Eight 
teams have provided a complete component separation pipeline 
capable to treat a realistic set of Planck temperature and polari- 
sation maps. Each methods differ in the external constraints they 
use, the physical modeling they assume and the algorithm they 
are based on. 

Among the available techniques, we choose to use the 
Generalized Morphological Component Analysis (hereafter 
GMCA) which is a blind method component separation method. 
In GMCA, each observation T° bs («) is assumed to be the linear 
combination of n c components {&'(n)}i = i t ...^ such that: 

T° bs (/z) = ^-s'W + n v («), (15) 

where n v (n) models instrumental noise. The general idea sub- 
tending this algorithm is the fact that the components, which 
result from completely different physical processes, have differ- 
ent spatial morphologies or structures. These morphological dif- 
ferences translate into a difference in their representation into 
a fixed waveform dictionary T>. If only a few coefficients of 
a fixed dictionary are enough to completely represent a given 
component, this component is said to be sparse in that dictio- 
nary T>. The dictionary succeeds in catching the general features 
that characterize the component morphology. That's why, sep- 
arating the observed map into components that maximize their 
sparsity in a given dictionary is an efficient strategy to distin- 
guish between physically different emission sources. In practice, 
a wavelet basis is a good choice for astrophysical components 
that overwhelmingly contain smooth spatial features. GMCA is 
a sparsity-maximization algorithm, a notion that we briefly in- 
troduce hereafter. Let {d-'(ra)} be the set of vector that forms the 
dictionary T>. Let aty = (s'(«), d J (n)) denote the scalar product 
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Fig. 2. Impact of the foregrounds residuals on the deflection field reconstruction on 12.5° x 12.5° square patches. Left: the input 
realisation of the deflection amplitude; Middle: the stack of 300 deflection estimates from the set I maps (synthetic Planck tem- 
perature maps - without any foregrounds residuals); Right: the stack of 300 deflection estimates from the set I-gmca maps (Planck 
temperature maps output of the GMCA component separation process). All maps have the same color table shown in the right-most 
part of the figure. 



coefficients between s'(n) and d ; (n). When D is an orthogonal 
wavelet basis, the following properties hold : 

<d , '(n),d''(H)> = Oif i + j 
<d'(ft),d''(»)> = 1 

s'(A) = 2<s , (n),d'(«))d-'(«) 



Then GMCA estimates the components {s'(«)} and the mixtures 
weights {Jl V i} by maxi mizing the sparsity of each component 
in D. As advocated in iBobin et all (120081) . a good sparsity es- 
timate is the sum of the absolute values of {a,y},j. Maximizing 
the sparsity of the components is then equivalent to minimizing 
this sparsity measure. The model parameters are estimated by 
GMCA as follows: 



mm 

{s'(a)},{A/l 



2|a i7 |s.t.||T° bs (n)-^[ vl s'(«)|| <e 



where e stands for the reconstruction error. The norm | 
for the usual (2 norm: 



T° bs (w) - yi VI s ! (n)|| = 



(16) 



stands 



(17) 



GMCA estimates the components {s'(ra)} which have only a few 
significant coefficients {07,} in the dictionary; i.e. the compo- 
n ents which are spar se in T>. Further technical details are given 
in lBobinet all (120081) . 

For Planck, the parameter e is chosen to be very small. In 
that case, the components {s'(«)} are estimated by applying the 
pseudo-inverse of the mixing matrix J[ to the observation chan- 
nels (T° bs («)}: 



s ! (n) = J^T° bs (/i) 



(18) 



where is the element at position {/, v) of the pseudo-inverse 

matrix of 31 defined as: J[ + = {jfjCj ' Jl T . Interestingly, the 
contribution of the component separation is then linear. As a 
consequence, the noise perturbing each component can be accu- 
rately known. Furthermore, the linearity of the separation guar- 
antees that the separation technique itself does not generate non- 
Gaussianity in estimated CMB map. Only the residual terms 



coming from the other components can create non-Gaussian fea- 
tures in the CMB. 

Another important consequence is that these properties give 
us a conservative method to estimate the point sources residu- 
als remaining within the CMB maps after foreground cleaning, 
as describe hereafter. Because each source has its own spectral 
property, component separation techniques fail at disentangling 
the point sources emission from the observed maps. As a re- 
sult, point sources remain mixed with the other components and 
the precise amount of the point sources emission by observation 
channels which has leaked in each component, is determined by 
the coefficients of the mixing matrix. More formally, in order to 
estimate the point sources residuals embedded in the foreground- 
cleaned CMB maps, quoted s ps («), one can apply Eq.[T8]to the 
simulated point sources in the observation channels {Ty S (n)}: 



s ps («) 



(19) 



where the elements {M^} form the column of the pseudo-inverse 
matrix which corresponds to the CMB component. Then the 
brighter point sources, which have been previously detected in 
the Planck-HFI channels, are masked out and the correspond- 
ing gaps are restored using an inpainting method. A detailed 
description of the mask and the restoration technique will be 
given hereafter in the Sect. |4] The final full-sky map we obtain 
is an estimate of the unresolved infrared point sources residuals, 
which contaminate the CMB temperature map after component 
separation with GMCA. This sky map is divided into patches of 
12.5 x 12.5 square degrees with 50% overlapping. We form a 
set of 300 point sources residuals square maps by selecting the 
patches with a maximal 30% masked area. 

We perform the component separation using GMCA on sets 
I-fg and Il-fg, each of 300 simulated patches generated following 
our idealized Planck sky model and described in the previous 
section (13. U . As an output of this process, we obtain two sets 
of 300 foreground-cleaned CMB temperature maps. Note that 
GMCA achieves the extraction of the foreground components as 
well. Unresolved point sources residuals are added to each map 
of these two sets. In the following, we will refer to the sets of 
lensed CMB maps with galactic dust, SZ effect and point sources 
residuals after the GMCA component separation, as sets I-gmca 
and II-gmca respectively. 
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3.3. CMB lensing reconstruction 

Here, we apply a discrete version of the quadratic estimator writ- 
ten by Okamoto & Hu (Eq. [6} on the different sets of simulated 
maps previously described, namely sets I, I-gmca, II and II-gmca. 
We seek to assess the foreground residuals impact on our capa- 
bility to reconstruct CMB lensing with Planck. 

3.3.1. Testing the estimator performances 

First, we give explicitly the expression of the discrete quadratic 
estimator that we derive from Eq. ©: 

hl(U) = ^P- V f(Hl)f(H 2 )^TT(«l,«2), (20) 

where U, u\ and «2 are wave-vectors related by u\ + «2 = U 
and A is the area of the sky patch considered. The sum is per- 
formed on the intersection of two disks T>\ and D^. The former 
is the zero-frequency centered disk defined by Ai|«i| < L max , 
where Aj denotes the frequency interval in the Fourier space ( i.e. 
the smallest nonzero positive frequency), which equals 2n/ VA. 
The latter, namely Di, is also a disk of radius L max but centered 
around U. It is therefore defined by Ai|(«2 - U)\ < L max . The 
wave-vectors «i, uo and U correspond to k\/A\, andL/Ai 
respectively. The normalisation Ajj(L) and the weighting func- 
tion Fjj(u 1,112) are respectively the discrete version of Eq. <j9j 
and Eq. (flOb . For the Planck-HFI experiment, we verify that the 
reconstructed potential field does not vary either we cut the sum 
in Eq. (1201 at L max = 2600 or we push it further. 

We study our capability to reconstruct a map of the inte- 
grated potential field with Planck-HFI idealized simulation, as- 
suming a perfect component separation without any foreground 
residuals. We apply the discrete quadratic estimator on set I 
maps (see Sect. 13.1b to obtain 300 estimates of the same real- 
isation of the projected potential field (p. Once stacking these 
estimates, the fin al map is an estimate of the input <p reali- 
sation. Following iHul (l2001bl) . we prefer to present our results 
in terms of the deflection field amplitude rather than the very 
smooth gravitational potential field, in order to highlight the in- 
termediate angular scales features. Fig. [2] shows the input de- 
flection field realisation, which has been used to simulate the 
lensing effect in the set I maps (on the first panel), as well as its 
reconstruction with the quadratic estimator applied on the set I 
maps (second panel). Even if the reconstruction noise is visible 
at smaller angular scales, the features of the deflection map are 
well recovered. 

Characterizing Planck sensitivity to the projected potential 
APS requires to account for both the CMB and the projected po- 
tential field cosmic variances. Thus, we move on to set II. As pre- 
viously, we apply the quadratic estimator (Eq. l20l i on the lensed 
CMB maps to reconstruct projected potential fields. Averaging 
over the APS of these individual <p field estimates gives an eval- 
uation of the quadratic estimator variance (as defined in Eq.[TT1i. 
The final reconstructed projected potential APS is obtained in 
subtracting the noise contributions, described in Sect. [2] from 
the variance. The former is related by Eq. (01 to the deflection 
APS shown in the Fig. [3] The error bars are estimated as the dis- 
persion between each individual deflection APS reconstruction. 
Thus, set II maps, which are idealized versions of the Planck- 
HFI sky assuming a perfect component separation, lead to a 
good reconstruction of the deflection APS up to L max = 2600. 
The error-bars evaluated here give an upper limit of the Planck- 
HFI sensitivity to the deflection APS. As one can see on Fig. |4] 



they are compatible with the theoretical Icr error-bars one can 
calculate from the Fisher formalism: 

ACf = + (21) 

where A^ff = 4n/LALA is a naive estimate of the independent 
available Fourier modes. The error bars estimated here will pro- 
vide us with a comparison level to quantify the impact of the 
foreground residuals. 

3.3.2. Impact of the foreground residuals 

Here we essentially redo the same analysis but using the full- 
simulation pipeline of our Planck-HFI demonstration model. 
The integrated potential field is extracted from the sets I-GMCA 
and II-GMCA described in Sect.|12l 

First, we aim at developing an intuition on the impact of fore- 
ground residuals on the deflection map reconstruction. We use 
the set I-gmca, in which both deflection field and foregrounds 
realisations are fixed. As previously, the 300 deflection field es- 
timates reconstructed with the quadratic estimator are stacked to 
produce an unique reconstructed deflection field shown in Fig. [2] 
We can see that the recovery of the underlying deflection field is 
still achieved even in presence of foreground emission and af- 
ter the GMCA analysis. The impact of the foreground residuals 
is nevertheless visible, mostly at angular scales larger than 2 de- 
grees whereas the intermediate angular scale features seem more 
preserved. 

For a more quantitative analysis, we move on to the impact 
of the foreground residuals on the deflection APS reconstruction. 
We use the set II-gmca (see Sects. [3TTI and l3.2b to ensure that the 
variances of the CMB, the deflection field and the foregrounds 
are accounted for. We reconstruct a deflection field estimate from 
each of the set II-gmca maps using the quadratic estimator given 
in Eq. d20t . Finally, we obtain the reconstructed deflection APS 
from the average variance over these 300 deflection APS esti- 
mates as described in Sect. 13.3.11 

The reconstructed binned deflection APS with the evaluated 
Icr errors is represented in the Fig. [3] Fig.|4]shows the difference 
between the reconstructed and the input deflection APS. First 
we report that the foreground residuals does not compromise 
the Planck-HFI capability to reconstruct the deflection APS - 
or equivalently the integrated potential APS. Fig. [3] shows that 
the APS reconstruction is preserved at the angular scales from 
L = 60 up to L = 2600. In this multipole range, the GMCA al- 
gorithm succeeds in letting unchanged the statistical properties 
of the lensed CMB temperature anisotropies, which suggests that 
this is a well-appropriated component separation tool for CMB 
lensing reconstruction. As for the first multipole bin, we report 
a 4<x excess of the deflection signal in the L — 2 to L = 60 
multipole range. We have checked that this bias is linked to the 
introduction of the unresolved point sources residuals in our sim- 
ulation pipeline. Interestingly, we find that this excess originates 
not from the level of residuals themselves, but mostly from the 
cutting procedural we use to extract the set of 300 square maps 
from our full-sky point sources residuals. We postpone a closer 
inspection of the low multipole lensing reconstruction behavior 



The cutting procedure involves a sphere-to-plan projection, an 
apodisation and a fit of the Fourier coefficients of the square map. From 
the tests we run (not presented here), the apodisation appears to have the 
most harmful impact on the high angular scale deflection reconstruc- 
tion. A complete study of the impact of the sphere to patches transition 
will be the subject of a companion paper. 
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Fig. 3. The deflection APS. Data-points are the binned APS 
reconstructed from Planck synthetic lensed CMB maps in 
two cases: (light blue/grey) the ideal case without any fore- 
grounds and (dark blue/black) the case with foregrounds residu- 
als from the GMCA output CMB maps. The fiducial deflection 
APS calculated with camb is figured by the (orange/solid) line; 
Orange/grey data points are the binned deflection APS estimates 
on the 300 input deflection field realisations. The horizontal and 
vertical intervals associated with the data points represent the 
averaging multipole bands and the lcr errors respectively. 

to the complete full-sky study (hereafter in Sect. [4j. Apart from 
the excess signal in the first bin, the impact of the foreground 
residuals on the deflection reconstruction is also slightly visi- 
ble at all angular scales, as seen in the Fig. [4] If the difference 
between reconstructed and input APS is still compatible with 
zero within the theoretical lcr errors in the 60 to 2600 multipole 
range, this residual bias appears more featured, more oscillat- 
ing than in the previous no foreground case. To quantify this 
degradation in the deflection APS reconstruction, we calculate 
the total error in unit of <x, defined as: 

A=Y|^1^|, (22) 

where C dA and C dd are respectively the reconstructed and in- 
put deflection APS in the b frequency band, and en, the lcr er- 
ror on C dd . Note that when evaluating the total error, we ex- 
clude the first bin bias, which has been discussed previously. 
With this definition, we find a total error, Aideai = 33, in the 
no foregrounds case whereas A emca = 41, in presence of GMCA 
foreground residuals. If the total error serves as quantifying the 
increase of the bias of the deflection APS reconstruction, one 
can also evaluate the increase of the errorbars. We find that 
the presence of foreground residuals results in a 10% increase 
of the errorbars on average. It suggests that at least an amount 
of the non-Gaussian foreground residuals is mixed up with the 
lensing signal in the reconstruction process. Foregrounds show- 
ing some small angular scale (around 15 to 5 arcmin) features - 
such as concentrated dust emission or brightest SZ clusters and 
point sources smeared by the instrument beam function - are po- 
tentially the more challenging for the CMB lensing reconstruc- 
tion. In addition, we remind that our sky model is intended to 
catch the dominant foreground features at the Planck-HFI fre- 
quencies. The sub-dominant existing emissions, as free-free or 



Fig. 4. Residual bias of the deflection APS reconstruction. Data- 
points figure the difference between the reconstructed and the 
input deflection APS averaging over 300 estimates: in the ideal 
case in light blue/grey and in the case with foreground residu- 
als in dark blue/black (consistently with the Fig. [3] caption). The 
lines show the non-Gaussian noise terms of the quadratic esti- 
mator, at the first-order (violet/dashed) and at the second-order 
(red/long dashed) in C Ad . Orange/grey colored band is the ana- 
lytical lcr error band derived from the Fisher formalism. 



synchrotron emissions, may marginally degrade our results as 
they are expected to slightly increase the foregrounds residuals 
whithin the CMB map. 

As a final remark, we note that, because of the first bin excess 
signal problem, linked to the sphere to patches transition, one 
might privilege a full-sky approach when seeking at a precise 
lensing reconstruction at the higher (L < 60) angular scales. 

4. Impact of masks on the full-sky lensing 
reconstruction 

4.1. Introduction 

From now, we move to a full-sky analysis of the CMB lensing 
effect. Some large area of the map, where the CMB signal is 
highly dominated by the foreground emission (e. g. the galac- 
tic plan, the point sources directions), have to be masked out. 
Cutting to zero introduces some mode-coupling within the CMB 
observables. As the lensing reconstruction methods rely on the 
off-diagonal terms of the CMB data covariance matrix, the map- 
masking yields some artifacts in the projected potential estimate 
if not accounted for. Several methods have been proposed to treat 
the masking effect for extracting the le nsing potential field from 
the WMAP data. In lSmith et all (I2007I) . the reconstruction is per- 
formed on the least-square estimate of the signal given the all- 
channels WMAP temperature data, requiring the inversion of the 
total data covariance matrix (S +N). However, as relying at least 
on the inversion of the noise covariance matrix, such an opti- 
mal data filtering approach is very CPU-consuming when ap- 
plied to the WMAP maps. The Planck-HFI resolution enables 
to provide 50 Mega pixels maps. Thus, the previous method to 
account for t he ma sking will be difficult to extend to Planck. In 
iHirata et aD, (l2008l) . the need for dealing with the noise covari- 
ance matrix is avoided by cross-correlating different frequency 
band maps. However, this method implies that no component 
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separation has been performed on the CMB maps before the 
lensing extraction. As a result, a lot of non-Gaussianities of fore- 
ground emission origin yields some artifact in the projected po- 
tential estimate, requiring a challenging post-processing to be 
corrected out. As previously mentioned, the Planck collabora- 
tion devotes strong efforts in the component separation activities 
and the c urrent developed m ethods have already proved their ef- 
ficiency dLeach et alj|2008l) . Moreover, the results we obtained 
with the demonstration analysis (see Sect. [3) tend to indicate 
that the lensing reconstruction is still doable after a component 
separation. Hence we plan to exploit the Planck frequency band 
maps to clean out the foreground emission before reconstruct- 
ing the lensing potential rather than using the cross-correlation 
based lensing estimator. As a consequence, we need an alter- 
native method to solve the masking issue in maps at the Planck 
resolution. Here, we propose to use an inpainting method, assess 
its impact on the CMB lensing retrieval and check its robustness 
to the presence of foregrounds residuals within the CMB map. 

First, we describe the hypothesis assumed and the tools we 
use to generate synthetic all-sky lensed temperature maps for 
Planck. Then we describe our full-sky lensing estimator and 
test its performances on some PLANCK-like temperature maps. 
Finally, we review the inpainting method and conclude in study- 
ing the effect of the inpainting on the projected potential APS 
reconstruction in two cases, first assuming a perfect component 
separation then in presence of point sources residuals. 

4.2. Full-sky simulation 

The formalism reviewed in Sect. [2] is almost fully applicable to 
the spherical case. In particular, the remapping equation (Eq. [TJ 
still holds, so that a lensed CMB sphere is given by: 



cision quality, LensPix is a well adapted tool for a CMB lensing 
analysis with Planck data alone. 

We obtained some Planck-HFI synthetic maps as the flat- 
sky case (see Sect. [3). White Gaussian noise realisations and 
Gaussian beam effect are added to the lensed temperature maps 
provided by the LensPix code. This Gaussian noise contribu- 
tion is fully defined by the all-channels beam-deconvolved APS 
given by: 



N 



1 



N, 



TT,v 



(24) 



T(h) = T(h + V^(m)), 



(23) 



where the V operator is to be understood as the covariant deriva 

tive on the sphere (iLewis & Challinorl l2b06). V(ft is identified to T^ hp \h) - V* — f/ m T/ m (n) 

be the deflection field d. Its zenithal and azimuthal coordinates 
can be calculated as the real and imaginary parts of a complex 
spin one field, using spin-weighted spherical h armonics t rans- 
form - the detailed calculation can be found in iHul (120001) and 
lLewii(l2005h 



where N? T ' V = (6 ivjhm o- T ) 2 exp 1)^/8^2], fwhm and 
<Xt are the full width at half maximum of the beam and the level 
of white noise per resolution element respectively, as given in 
Table Q] The noise map is generated using the map creation tool 
of the HEALPix package. 

For the lensing reconstruction analysis, we prepare two sets 
of 50 all-sky maps with 1.7 arcmin of angular resolution (the 
HEALPix resolution parameter nside = 2048). In each set, maps 
are the lensed CMB temperature plus the Planck-HFI nominal 
Gaussian noise. 

4.3. Full sky lensing reconstruction 

We carry out an integrated potential estimation tool based on th e 
full-sky version of the quadratic estimator derive d inlHul (1200 1 ah . 
We cl osely follow the prescription given in lOkamoto & Hul 
(120031) to build an efficient estimator, so that: 

n (0) r 

4>lm = L + - J dh (T (h ''\h)VT <H \h)) ■ VY* LM (fi), (25) 

where and T (w) are respectively high-pass filtered and 

weighted lensed CMB temperature field, given by: 



C 1 

lm 

cl 



(26) 



r (w) («) = £ _L_f te y te (a) 



The LensPi^0 package described in ILewis! d2005l) aims at 
generating a set of lensed CMB temperature and polarisation 
maps from the analytical auto- and cross-APS of {T, E, B,<p}, 
the temperature, the E and B polarisation modes and the line- 
of-sight projected potential respectivel y. The maps are pro vided 
in the HEALPixQ pixelisation scheme dGorski et al.ll2005l) . The 
CMB lensing simulation is achieved in remapping the anisotropy 
fields according to Eq. ( f23b of a higher resolution map using a bi- 
cubic interpolation scheme in equi-cylindrical pixels. A lensed 
temperature map at the Planck resolution (nside = 2048), can 
be computed in about 5 minutes on a 4-processors machine. The 
relative difference between the lensed temperature APS recon- 
structed on su ch a map and the a nalytical lensed APS obtained 
with CAMBQ (ILewis et al.ll2000T) is below 1% up to / = 2750. 
Here, to conservatively ensure a relative error below 1%, we 
choose a multipole cut at L max = 2600. By its speed and pre- 



The covariant derivative operator V applied on the spherical har- 
monics can be expressed in term of the spin +1 projectors e ± - 
eg ± ie<t, and the spin weig hted spherical harmonics, as explained 
in lOkamoto & Hul (120031) . The quantities appearing in Eq. (f25b 
can thus be calculated by direct and inverse (spin zero) spheri- 
cal harmonics and spin-weighted (spin +1) spherical harmonics 
transforms. As for the estimator normalisation, quoted A^°\ it 
identifies with the Gaussian contribution to th e estimator vari- 
ance; F or its expression, we refer to Eq. (34) in lOkamoto & Hul 
(120031) . Then, extending to the spherical case the calculations re- 
viewed in Sect. [2] the covariance of the integrated potential field 
estimator (piM, averaged over an ensemble of CMB and gravita- 
tional potential fields realisations, depends on the potential APS, 
so that: 



('Plm'Pl'M'} - &LL'<>MM' 



(c* 



(27) 



5 http : //cosmologist . info/lenspix/ 

6 The acronym for 'Hierarchical Equal Area isoLatitude Pixelization' 
of a sphere (see http: //healpix. jpl .nasa. gov/index . shtml ). 

7 The 'Code for Anisotropics in the Microwave Background' is a so- 
called Boltzmann's code described at http : //camb . info/ 



where N®\ N^' and are zeroth, first and second order in 
Cjf noise terms respectively. To calculate the first sub-dominant 
noise term, one c an use the expressi on derived in the flat-sky 
approximation bv lKesden et al.l(l2003l). L ikewise, the second or- 
der term is given in lHanson et al.l (120091) . As these terms are an 



r(l) 



,<2> 
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Fig. 5. All-sky maps in Galactic Mollweide projection. Upper panels: (left) a Planck-HFI synthetic CMB temperature map in 
milliKelvin, (right) the union mask defined in Sect. 14.51 The grey region shows the observed pixels whereas in black are the rejected 
ones. Lower panels: left panel shows the restored CMB map obtained by applying the inpainting process (described in Sect. l4.4b on 
the previous upper-left CMB map masked according to the union mask and the right panel shows the difference between the input 
(upper-left) CMB map and the restored (lower-left) one. 



order of magnitude smaller than the dominant A^ 0) term and be- 
cause the flat-sky approximation is known to be robust fo r the 
potential estimator noise calculation dOkamoto & Hull2003l) . we 
assume and will verify that the deviation due to the flat-sky ap- 
proximation is negligible. 

Within the framework of Monte-Carlo analysis, we build a 
projected potential APS estimator so that: 



C 



1 N \ 1 — 

2~ltiX | ^ m|2 

i=i M 



(28) 



where 4>' LM is the integrated potential field estimate on the i 
CMB temperature realisation, ; e {1, . . . A'}. Note that since 

(2) 

and iV£ depend on the potential APS itself, it should be eval- 
uated and subtracted iteratively. Here we calculate it once from 
the theoretical integrated potential APS. 

Finally, we test our APS estimator on a set of 10 lensed CMB 
temperature maps of 50 millions of pixels, including the nomi- 
nal Planck noise, generated as described in Sect. 14.21 As in the 
flat-sky case, sums in the spherical harmonic space are cut at 
L max =2600. The results, compiled in the form of an integrated 
potential APS estimate averaged over the 10 trials (see Eq. d28b), 
are shown in the left panels of Fig. [6] 

4.4. Inpainting the mask 

We choose to take into account the cutting effect of the temper- 
ature map before any lensing reconstruction rather than making 
any changes in the quadratic estimator (given in Eq. [25] ) to ac- 
count for the mask. This approach is motivated by the fact that 
the high quality and the large frequency coverage of the Planck 



data allow to reconstruct the CMB temperature map on roughly 
90% of the sky. It therefore suggests that a method intended to 
fill the gap in the map can be applicable. 

Several of such methods, referred to as inpainting, 
have been recent l y dev eloped since the pioneering work of 
iMasnou & Morell (fl998h . The general purpose of these meth- 
ods is to restore missing or damaged regions of an image to 
retrieve as far as possible the original image. For the CMB lens- 
ing reconstruction, the ideal inpainting method would lead a re- 
stored map having the same statistical properties than the un- 
derlying unmasked map. To use a notion briefly mentioned in 
the section 13.21 the masking effect can be thought of as a loss 
of sparcity in the map representation: the information required 
to define the map has been spread across the spherical harmon- 
ics basis. That's why, the inpainting process can also be thought 
of as a restoration of the CMB temperat ure field sparcity in a 
conveniently chosen waveform dictionary. lElad et al] d2005l) in- 
troduced a sparsity-based technique to fill in the missing pix- 
els. This me thod has been extended to the sphere in lAbrial et al.l 
(I2008L [2007D . In a nutshell, the masked CMB map is modeled as 
follows : 



T{h) = M(h)T(h) 



(29) 



where Ai(h) stands for a binary mask the entries of which are 
one when the pixel is ob s erved and zero when it is missing. As 
emphasized in Elad et al.l d2005l) : lAbrial et al.l d2008l) . if T(h) has 
a sparse representation in a given waveform dictionary T> (see 
13.2b . masking is likely to degrade the sparsity the CMB map in 
T>. Let {d'{ft)\ be the set of vector that forms the dictionary T>. 
Let aj denote the scalar product (so-called coefficients) between 
T(h) and d-i(h) : aj = (T(h), c/ ; (n)). For the sake of simplicity, 
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Fig. 6. Impact of the masking corrected by an inpainting process. The left panels are for the full-sky Planck synthetic lensed 
temperature maps whereas the right panels compile the results drawn from the masked lensed temperature maps restored with 
the inpainting process described in Sect. I4.4I Upper panels: the reconstructed deflection APS C dd , obtained from the projected 
potential APS estimate of Eq. d28l ). The (light blue/grey) crosses show the C dd per multipole and the (dark blue/black) data points 
are the band-power C dd . The horizontal intervals represent the averaging multipole bands and the vertical ones the lcr error. The 
(orange solid) line figures the fiducial deflection angle APS C dd and the (green dashed) line the total noise of the quadratic estimator 
(A^ 0) + A^ 1 ' + A^ 2) ). Lower panels: the bias of the deflection APS reconstruction AC dd defined as the difference between the 
reconstructed deflection APS C dd and the input deflection APS C dd . Consistently with the upper panels, (light blue/grey) crosses 
are AC dd per multipole and (dark blue/black) data points are the band power AC dd with the associated averaging multipole widths 
(horizontal intervals) and lcr error (vertical intervals). The large (orange/grey) band shows the analytical +l<x errors per band 
power expected for the quadratic estimator (see Eq. l2Tb. Finally, the lines show the non-Gaussian sub-dominant noise terms at the 
first-order (violet/dashed) and at the second-order (red/long dashed) in Cfj . 



we further assume the set \d^(h)} forms an orthonormal basis. 
Recovering the missing pixel can then be made by for a solu- 
tion that minimizes the sparsity of the T(h) in G. As in 13.21 an 
appropriate sparsity estimate of T(h) in T> consists in measur- 
ing the sum of the absolute values of a, = (T(h), d J '(h)). The 
recovered CMB map is then obtained by solving the following 
optimization problem : 

min V \(T(h),d j (h))\ s.t. \\T(h) - M(h)T(h)\\ < e (30) 

{T(h)\ 

J 

w here e stands for th e reconstruction error. It has been shown 
in lAbrial et alj d2008l) that this inpainting technique leads to very 
good C MB recovery r esults. Our i npainting algorithm can be 
found in lAbrial etall d2007l 120081) and its implementation is 
based on the Multi-Resolution on the Sphere (MRS) packag^l 

8 http://jstarck.free.fr/mrs.html 



In the following, we seek to assess the impact of the inpainting 
mask correction in PLANCK-like maps on a CMB lensing recon- 
struction. 

4.5. Effect of inpainting 

First, we have to choose a realistic mask, which could apply to 
the forthcoming Planck temperature map. However, depending 
on the details of the component separation pipeline the different 
metho ds developed in the Planck consortium (see iLeach et"aD 
120081) . yield to slightly different masks. Furthermore, the mask 
size is not a necessary criteria for the final choice of the compo- 
nent separation method that will be selected for the Planck data 
analysis. We thus adopt a conservative approach, which consists 
in choosing the union of the masks provided by each of the meth- 
ods at the time of the Component Separation Planck Working 
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Group second challenge dLeach et al.ll2.Q08h . Such a mask, here- 
after referred to as the union mask, rejects about 1 1% of the sky, 
as shown in the Fig [5] 

Then the 10 PLANCK-like lensed CMB temperature maps we 
have generated (see Sect. 14.2b are masked according to the union 
mask and then restored by app lying the inpainting method de- 
scribed in Abrial et alj d2008l) . From each of these mask cor- 
rected maps, we extract a projected potential field using the 
quadratic estimator of Eq. d25l l. As previously, the results are 
compiled in the form of the average projected potential APS Cfj 
following Eq. (128b . The reconstructed deflection APS, given by 
C dd = L(L + 1)C£ , as well as the bias between estimated and 
fiducial deflection APS AC dd are shown in the right panels of 

Fig. m 

We find that the mask corrected by the inpainting results 
in a marginal increase (~ 4%) of the lcr errors on the esti- 
mated deflection APS (hence on the projected potential APS). 
Masking and inpainting causes an increase of the reconstructed 
APS bias AC dd arising mostly at large angular scale corre- 
sponding to multipole L < 300. However, this bias is weaker 
than the sub-dominant second-order in cfj non-Gaussian bias. 
Fig. [6] shows a clear increase of power in the very first multi- 
pole band (2 < I < 10). In this multipole range, Planck is 
not e xpected to achieve a g ood reconstruction of the potential 
APS dHu & Okamotoll2002l . From the multipole L = 300 up 
to L — 2600, the bias stays below the first-order in Cy non- 
Gaussian bias and is compatible with the theoretical lcr errors 
expected for the quadratic estimator. From fully controlling the 
inpainting impact, one might want to push further the study by 
analytically calculating or Monte-Carlo estimating the mask in- 
duced bias. However, it is not mandatory for reconstructing the 
projected potential APS with Planck. The masking effect, once 
corrected by inpainting, becomes a sub-dominant systematic ef- 
fect that can be safety neglected. 

4.6. Robustness against the unresolved point sources 

Up to now, we have handled independently two important issues 
linked to the presence of foreground emissions in the observa- 
tion maps, the impact of the foreground residuals after compo- 
nent separation with GMCA in Sect. 13.31 and the impact of the 
masking corrected with the inpainting method in the previous 
subsection (Sect. |4~5l l. We found that none of them compromises 
our ability to reconstruct the deflection APS. In a more realistic 
approach, these two issues should be handled altogether, as the 
inpainting process is intended to be applied on a CMB map con- 
taminated by foreground residuals. The presence of foreground 
residuals is susceptible to harden the inpainting process and con- 
sequently degrade the CMB lensing recovery. Here, we assess 
the robustness of the deflection reconstruction on masked and 
inpainted CMB maps when adding infra-red point sources resid- 
uals. This choice is motivated for two reasons. The point sources 
residuals after component separation is a well-known matter of 
concern in any CMB non-gaussianities analysis and the emis- 
sion of the infra-red sources population is one of the major fore- 
ground contaminant at the Planck-HFI observation channels. 

We use the full-sky map of infra-red point sources resid- 
uals after a component separation using GMCA, we had esti- 
mated in Sect. 13.21 This point sources residuals is added to the 
10 synthetic Planck lensed CMB temperature maps described 
in Sect. 14.21 Then we redo the same analysis than previously 
in Sect. 14.51 the union mask is applied to the maps, cutting out 
the brightest infra-red sources, which have been detected during 



the Co mponent Separatio n Planck Working Group second chal- 
lenge dLeach et al.ll20 08). The 10 masked maps are restored us- 
ing the inpainting method before being ingested in the full-sky 
quadratic estimator of the projected potential field. The results 
of the whole analysis are presented in the form of the average 
reconstructed deflection APS and bias, and shown in Fig. [7] 

We find that the inpainting performances are only marginally 
degraded (at < lcr level) by the presence of point sources resid- 
uals within the CMB maps, and this degradation occurs mainly 
at the two multipole extremes. At the lower multipoles (L < 30), 
the APS deflection reconstruction suffers from a lcr increase of 
the bias, whereas at the higher multipoles, only the error bars 
increase. We conclude that the inpainting method succeeds in 
keeping the statistical properties of the CMB map unchanged 
even in presence of highly non-Gaussian foreground residuals 
and it is a qualified method to handle the masking issue when 
seeking at a CMB lensing recovery. In addition, the results com- 
piled in the Fig. [7] give the total impact of point sources on the 
deflection reconstruction, as they account for both the masking 
of the bright detected sources and the unresolved residuals. We 
report that point sources are responsible for a total 13% increase 
of the lcr errors on the reconstructed APS deflection, mainly 
induced by the unresolved residuals. As a summary, the major 
nuisance of point sources is related to the masking of the bright 
ones, which tend to increase the bias on the reconstructed deflec- 
tion APS, whereas the unresolved residuals results mainly in an 
increase of the errors on the deflection retrieval. 

Conclusions 

The High Frequency Instrument (HFI) of the Planck satellite, 
which has been launched on the 14 th of May 2009, has the sensi- 
tivity and the angular resolution required to allow a reconstruc- 
tion of the CMB lensing using the temperature anisotropics map 
alone. The pioneer works to put evidence of the CMB lensing 
within the WMAP data are not directly applicable or not well- 
optimized to the Planck data. First, one might want to take ben- 
efit of the efficient component separation algorithms developed 
for Planck before applying a CMB lensing estimator rather than 
to correct the lensing reconstruction from the bias due to the 
foreground emission afterward. Second, we need an efficient and 
manageable method to take into account the sky cutting within 
the 50 Mega-pixels maps provided by Planck. In addition, the 
CMB lensing is related to another burning thematic: character- 
izing the non-Gaussianities of the temperature anisotropies (pri- 
mordial non-Gaussianities, cosmic string, etc.). 

We have implemented both the flat-sky and the all-sky ver- 
sions of the quadratic es t imator of the projected potential field 
described in iHul d2001bb : lOkamoto & Hul j2003h to apply them 
on Planck synthetic temperature maps. First, within the flat-sky 
approximation, we have prepared a demonstration model, which 
consists in r unning GMCA, a component separation method 
described in iBobin etaTI d2008l) . on Planck frequency chan- 
nel synthetics maps, containing the lensed CMB temperature, 
the Planck nominal instrumental effects (modeled by a white 
Gaussian noise and a Gaussian beam) and the three dominant 
foreground emissions at the Planck-HFI observation frequen- 
cies, namely the SZ effect, the galactic dust and the infra-red 
point sources. We have performed a Monte-Carlo analysis to 
quantify the impact of the foreground residuals after the GMCA 
on the projected potential field and APS reconstructions. Then, 
we ha ve moved on to the full-sky case, using the LensPix algo- 
rithm dLewisl2005l) to generate lensed CMB temperature maps at 
the Planck resolution. We have performed a Monte-Carlo anal- 
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Fig. 7. Robustness of the inpainting to the unresolved point 
sources. The upper panel shows the reconstructed deflection 
APS whereas the lower panel the bias on the deflection APS re- 
construction. The (dark blue/black) data points show the result of 
the full-sky quadratic estimation on the set of lensed CMB maps 
with point sources residuals, which have been masked then re- 
stored with the inpainting technique, as described in Sect. I4.6I 
For comparison, the results obtained in Sect. I4.51 in the case 
without any sources residuals, have been copied out on here as 
(light blue/grey) crosses. The fiducial analytical deflection APS 
as well as the total reconstruction noise and the two APS bias 
are shown following the same representation code as in Fig. |6l 
and horizontal and vertical intervals have the same meaning as 
described in the Fig.[6]caption. For a sake of readibility, only the 
band power reconstructed APS and APS bias are represented, 
and the theoretical error bars by multipole bins are not shown. 



ysis to tackle the m asking issue; we ha ve used the inpainting 
method described in lAbrial et~ai] d2008l) to restore the Planck 
synthetic temperature maps, masked according to a realistic cut 
out of 1 1 % of the sky, accounting for the bright detected point 
sources. By applying the projected potential quadratic estimator 
on these restored maps, we have studied the impact of the in- 
painting of the mask on the Planck sensitivity to the projected 
potential APS. Finally, we have assessed the total impact of the 
point sources emission, in confronting the inpainting method to 
the unresolved point sources residuals. 

Results 



1 . Within our flat-sky demonstration model, we found that the 
reconstruction of the projected potential field is still feasible 
after a component separation using GMCA. More quantita- 
tively, the foreground residuals in the GMCA output CMB 
maps lead to a 10% increase of the lcr errors on the projected 
potential APS reconstruction when applying the quadratic 
estimator. The GMCA process results in an increase of the 
dispersion of the projected potential APS reconstruction, but 
this dispersion remains within the theoretical lcr errors at 
all angular scales but the L < 60 multipoles, in which the 
flat-sky analysis is expected to show some limitations any- 
way. Such a study dealing with the impact of a component 
separation process on the CMB lensing reconstruction had 
never been performed before. Our results allow us to assess 
that applying a component separation algorithm on the fre- 
quency channel CMB maps before any lensing estimation is 
a well-adapted strategy for the sake of the projected potential 
reconstruction within Planck. 

2. For the full-sky reconstruction of the projected potential APS 
with Planck, we report that a realistic 1 1% of the sky mask, 
applied on some PLANCK-nominal lensed CMB temperature 
maps, has a negligible impact on the CMB lensing signal 
retrieval process, w henever it has been corrected by the in- 
painting method of lAbrial et alj d2008l) beforehand. More 
precisely, the bias on the estimated projected potential APS 
induced by the mask after inpainting is always either com- 
patible with the theoretical lcr errors (from I = 300 up to 
I = 2600) or weaker than the second-order in cf^ non- 
Gaussian bias (in the / < 300 range). The major impact 
of the inpainting correction on the projected potential APS 
arises at the larger angular scales (2 < I < 10), which are not 
expected to be well-reconstructed with Planck. In addition, 
these results have not significantly changed after the intro- 
duction of unresolved point sources residuals. When treating 
the point sources emission in a comprehensive way, we re- 
port a 13% increase of the lcr errors on the reconstructed 
deflection APS on average, resulting mainly from the un- 
resolved point sources residuals, whereas the level of bias 
is marginally increased at low multi poles. We conclude that 
applying the inpainting method of lAbrial et alj d2008l) be- 
forehand is a good strategy to take into account the masking 
issue when seeking at reconstructing the projected potential 
with Planck. 



Perspectives Our results on the CMB lensing reconstruction 
are the first step to elaborate a complete analysis chain dedicated 
to the projected potential APS reconstruction with the Planck 
data. Such a CMB lensing reconstruction pipeline should involve 
a component separation and a bright point sources detection fol- 
lowed by an algorithm to correct from the mask (e.g. the in- 
painting method) before applying a quadratic estimator of the 
projected potential field on the resulting CMB temperature map. 

We plan to go on developing in parallel both the flat-sky and 
the full-sky reconstruction tools. The flat-sky tools will allow us 
to perform a multi-patches CMB lensing reconstruction in cut- 
ting several hundred patches out of the most foreground cleaned 
region of the full-sky map. Using such a method requires a quan- 
titative study of the impact of the sphere-to-plan projection and 
the sharp edge cuts effects beforehand. As a first task, we should 
test whether our results concerning the feasibility of reconstruct- 
ing the CMB lensing after a component separation and after an 
inpainting of the mask still hold when dealing with the fully real- 
istic Planck simulation (including e.g. non axisymmetric beam, 
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inhomogeneous noise and correlated foreground emissions). As 
long as we can demonstrate we have a sufficient control on sys- 
tematics, we will be ready to measure the projected potential 
APS with Planck alone. Such an additional cosmological ob- 
servable is expected to enlarge the investigation field accessible 
to the Planck mission from the primordial Universe to us. 
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